LSND anomaly from CPT violation in four-neutrino models 
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The LSND signal for — > Ve, oscillations has prompted supposition that there may be a fourth 
light neutrino or that CPT is violated. Neither explanation provides a good fit to all existing 
neutrino data. We examine the even more speculative possibility that a four-neutrino model with 
CPT violation can explain the LSND effect and remain consistent with all other data. We find that 
models with a 3 + 1 mass structure in the neutrino sector are viable; a 2 + 2 structure is permitted 
only in the antineutrino sector. 



I. INTRODUCTION 

The LSND experiment has found evidence for v^^ —>■ 
oscillations at the 3.3(7 level 0,13, with indications for 
— > oscillations at lesser significance 0, 01- The 
combination of the LSND data with the compelling ev- 
idence for oscillations in solar, atmospheric, accelerator, 
and reactor neutrino experiments cannot be adequately 
explained in the standard three-neutrino picture with 
CPT conservation 4]. Extensions to models with four 
light neutrinos (with the extra neutrino being sterile) |4j 
or CPT violation 0, 1^ with three neutrinos have been 
proposed to accommodate all neutrino data. However, in 
both cases, recent analyses indicate that neither scenario 
provides a good description of the data j3, 13 . 

The MiniBooNE experiment ^ is now taking data that 
will test the LSND oscillation parameters^ in the 
Ue channel. A positive result in MiniBooNE will rule 
out current versions of CPT-violating models, while a 
negative result will rule out four-neutrino models with 
CPT conservation. In either case, the surviving models 
will still not give a good fit to all data. 

In this letter we consider the very speculative possibil- 
ity of CPT violation in Jour-neutrino models. The CPT 
violation is manifested as different mass matrices for neu- 
trinos and antineutrinos^ . We find that such a scenario 
can accommodate the data only if the sterile neutrino is 
weakly coupled to active neutrinos. Thus, while 3 -I- I 
models are viable, 2 + 2 models (in which the sterile 
neutrino is strongly coupled to active neutrinos in solar 
and/or atmospheric oscillations), are not. A hybrid so- 
lution (3 + 1 for neutrinos and 2 + 2 for antineutrinos) 



^ The bulk of the parameter region allowed by — > data 
from LSND and KARMEN Op and Pe data from 

the Bugey reactor experiment lies in a narrow band in 

(sin^ 2eL,5m\) space along the line described approximately by 
sin^ 2ei((5ni2 )i-64 = 0.0025 between 5m\ ~ 0.2 and 1 eV^. A 



small allo wed region near Sm?^ ~ 7 eV'^ and siv?29]^ = 0.004 
also exists Il2i . 

Whether such a model can be constructed using nonlocality of 
the interactions with out violating Lorentz invariance is still a 
matter of debate Il3l . 



is also possible. 



II. FOUR NEUTRINOS OR CPT VIOLATION? 

There are two types of four-neutrino models: (a) 3 + f , 
where active neutrinos have mass-squared differences and 
mixings similar to the standard three-neutrino model 
that describes solar and atmospheric data, and (b) 2 + 2, 
where there are two pairs of closely spaced mass eigen- 
states, one of which accounts for the solar neutrino data 
and the other for the atmospheric neutrino data. The 
3+1 models with CPT conservation are disfavored be- 
cause the Bugey reactor and CDHSW accelerator Q 
experiments put constraints on oscillation amplitudes for 



and 



survival, respectively, which to- 



gether imply an upper limit on the LSND oscillation am- 
plitude that is below the experimental value 0, Hsl ll^ . 
The 2 + 2 models with CPT conservation are ruled out 
because the combination of solar and atmospheric data 
do not allow enough room for a full sterile neutrino 
To resurrect the scenario in which sterile and active states 
are only weakly coupled, it has been shown that extend- 
ing 3 + 1 models by an extra sterile neutrino improves 
the fit to short-baseline data substantially 0| . The con- 
straints on 2 + 2 models may be relaxed by including 
certain small neglected mixing angles in the analysis ^I8| . 

In CPT-violating models with three neutrinos, the 
mass spectra, and hence the mass-squared differences, 
are different for neutrinos and antineutrinos'^. In the 
original versions of these models , the neutrino mass- 
squared differences accounted for the solar and atmo- 



spheric oscillations (but not the weak v^^ 



signal 



in LSND), while the antineutrino mass-squared differ- 
ences accounted for the LSND and atmospheric oscilla- 
tions. With the addition of the KamLAND data indicat- 
ing Ve —> i>e oscillations at the solar Sm? scale, the an- 
tineutrino mass-squared differences were adjusted to ac- 
count for the LSND and KamLAND oscillations (but not 



^ A comparison of s olar and reactor neutrino data can constrain 
CPT violation E3; a 90% C. L. limit of |<5m2-5m2| < 1.3 xlO-^ 
eV^ was found in Ref. 1201. 
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oscillations of atmospheric antineutrinos) An anal- 
ysis of the atmospheric data showed that the modified 
CPT-violating model did not give a good fit and is ex- 
cluded at the 3(7 level Q. Thus, neither four- neutrino 
models with CPT conservation nor three-neutrino mod- 
els with CPT violation provides a consistent explanation 
of all the data including LSND. 



III. FOUR NEUTRINOS WITH CPT 
VIOLATION 

The neutrino flavor states Va {a = e, /z, r, s) are related 
to the mass eigenstates Vi by a unitary matrix U , with 

4 

i=l 

If CPT is not conserved, then the corresponding unitary 
matrix U for antineutrinos, given by 

P„ = E IJaiVi , (2) 
i 

will not necessarily be equal to U. Furthermore, the 
neutrino and antineutrino eigenmasses will not neces- 
sarily be the same. We will assume that there are 
three different mass-squared difference scales for neu- 
trinos, (Sm^ ^ Sm^ ^ S'tti'l, that can explain the so- 
lar, atmospheric and LSND data, respectively, and that 
the corresponding mass-squared differences for antineu- 
trinos are similar to those for neutrinos, i.e., Sfn^ « <5m^, 
Sfri^ « Snr^, and Srfij^ « Smj^. 



A. 3 + 1 models 

In 3 -I- 1 models there is one neutrino mass well- 
separated from the others by <5m|^, and the sterile neu- 
trino couples strongly only to the isolated state. There 
are four different mass spectra in 3 -I- 1 models, depend- 
ing on whether the isolated state is above or below the 
others, and whether the other three neutrino states have 
a normal or inverted mass hierarchy. We consider the 
case with > mi,TO2,TO3 and normal hierarchy, which 
implies Smli ~ (5m|2 — Srnl^ — 5m\ ^ 5m\i ~ 5m\2 — 
Sm^ ^ 5m2i = Srn^, with similar relations for Sfrifj (the 
argument is similar for the other cases). Then the rel- 
evant oscillation probabilities for neutrinos are approxi- 
mately, 

Pi^^e -> J^e)solar ^ 1 " 4|C/ei |[/,2 sin^ A, , (3) 

^ i'/^)atm ^ l-4|[/^3|'(l-|C/M3l'-|C^M4nsin2 A,, 

(4) 

P(l/^ -» I^e)MiniBooNE ^ ■i\Ue4\'^\U^4:\'^ sin^ Al , (5) 



^ f^M)cDHSw - 1 - 4|i7^4|'(l - \U^4\^) sin^ Al , 

(6) 

where A^; = Srn^L/ {AEi,) is the usual oscillation argu- 
ment for X — a, or L. The relevant oscillation proba- 
bilities for antineutrinos are approximately, 

P{Ve ^ j>e)KamLAND ^ 1 - 4|i7ei ^ |C7e2 P Sin^ A, , (7) 
P{l>^. ^ t>;.)atm ^ 1 " 4| C7^3 P ( 1 " I I ' " I H siu^ A,, 

(8) 

'^e)LSND ^ 4|C7e4r|i7;.4p siu^ Ai , (9) 
Pi^e ^ l>e)Bugoy ^ 1 " 4| C7e4 ^ (1 - | C/e4 P) siu^ Ai , (10) 

where A^, = 5mlL/{-iEi,). We note that LSND has al- 
ready made a measurement for PminiBooNE, although at 
the 2(7 level it is consistent with both Plsnd and zero. 

If CPT is not conserved^, then in general U ^ U . 
The Bugey-I-CDHSW bound on LSND_and MiniBooNE 
can now be evaded since Bugey limits | Uei \ and CDHSW 
limits |C/p4|, but \U^4\ and |J7e4| are no longer tightly 
constrained. The bounds on the amplitudes for — > Dg, 
and — > Ve oscillations must be determined separately^: 

(i) The best constraint on \Ue4\ comes from Vf, disap- 
pearance performed during GALLEX testing with 
a ^^Cr neutrino source [23 ; it does not extend as 
low in 5m\ as the Bugey constraint and is about 
an order of magnitude weaker at high 5m\. (The 
GALLEX survival probability is given by Eq. [Hil 
with Uei replaced by Uei). When combined with 
the limit on |C/^4| from CDHSW 14|, a range 
of oscillation amplitudes (4|C/e4p|J7^4p) 

is excluded; see Fig. A positive MiniBooNE 
result that exceeds the upper bound in Fig. 
will rule out 3-1-1 models with CPT violation, 
while a null or positive result that obeys the up- 
per bound is easily accommodated. A MiniBooNE 
measurement consistent with the region allowed by 



* The reason that 3 + 1 models with CPT conservation (where 
U = U) are disfavored is that |(7^4p mu st be small from a 
combination of the strict CDHSW limit |l4l and the large mix- 
ing of atmospheric u^^, and |C/e4p must be small from a com- 
bination of the strict Bugey limit llll and the large mixing 
of solar Ve and KamLAND 9^, lea ding to an upper bound on 
the LSND amplitude A\UeA\'^\U^A? [11111. The CDHSW and 
Bugey limits are i5r?i|^ dependent (as is the LSND allowed ampli- 
tude), but a comparison at all i5rri|^ shows that nowhere does 
the LSND+KARMEN 95% C. L. allowed region overlap the 
95% C. L. allowed region from the other experiments |21| . The 
most recent comprehensive analysis concludes that the 3 + 1 
models have a goodness of fit of at most 5.6 X 10~^ 01 ■ 

^ Although neutrino telescopes are capable of probing the LSND 
scale 1231 . they are unable to test CPT-violating schemes be- 
cause neutrino and antineutrino oscillation probabilities are not 
measured separately. 
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FIG. 1: Upper bound (solid) on the i^^ — » oscillation am- 
plitude 4|f/e4p|C/M4|^ from the GALLEX limit on \Ue^\ and 
the CDHSW limit on \Uy,^\ (90% G. L. results are used in 
both cases). The dot-dashed line is the 99% G. L. u ppe r 
bound from Bugey and GDHSW if CPT is conserved |21| . 
Also shown are the expected sensitivity (dashed) of the Mini- 
BooNE experiment and, for comparison, the allowed region 
(within the dotted lines) for 4jf7e4p|f/M4p from a combined 
analysis of LSND and KARMEN data, both at the 90% 
G. L [13. 

LSND and KARMEN would mean that the values 
of 4|;7e4p|?/^4|^ and 4|C7e4p|t//i4p are similar, but 
the Bugey-I-CDHSW bound demands IJ ^ U . In 
fact, a MiniBooNE measurement that lies below 
the GALLEX-I-CDHSW upper limit and above the 
Bugey-|-CDHSW limit, would place a lower bound 
on the amount of CPT violation. 

(ii) The best accelerator constraint on |C/^4| comes from 
the CCFR experiment that searched for j/^ disap- 
pearance |2^]; there is no limit for (5m| < 7 eV^. 
Although I C7^4 1 cannot be so large as to disturb the 
usual fits to atmospheric data (which indicate the 
dominant oscillation is at the Jm^ scale), 1^7/^41^ 
can probably be of order a few per cent, which al- 
lows the LSND amplitude 4|f7e4p|C/;i4p to be large 
enough to account for the LSND v^^ data, at 

least for the larger values of 8m\^ allowed by LSND 
(near 1 eV^ and 7 eV^); see Fig.|21 

We note that fits for neutrinos and antineutrinos no 
longer must agree, so that, e.g.-, the solar fits are now in- 
dependent of the KamLAND fits. Thus, any differences 
that might occur between the KamLAND and solar neu- 
trino allowed regions could be explained by CPT vio- 
lation. Similarly, results for atmospheric neutrinos and 
antineutrinos (which could perhaps be measured sepa- 



FIG. 2: Lower bounds on 4|0'fi4p(l - |'^m4|^) (the amplitude 
for atmospheric i/^ survival at the LSND mass scale) from 
the Bugey limit on Ve disappearance and the t-^ — > Ve oscilla- 
tion amplitude indicated by LSND and KARMEN (90% G. L. 
results are used in both cases). 

rately in MINOS could also be different and stiU be 
easily accommodated in the model. On the other hand, 
if there were no discernible difference between the neu- 
trino and antineutrino fits to solar, KamLAND, and at- 
mospheric data, the amount of CPT violation needed 
would be small, since then the only differences between 
neutrino and antineutrino parameters would occur in the 
small mixings between the sterile and the active states. 



B. 2 + 2 models 

Constraints on 2 -I- 2 models are different from those on 
the 3 -t- 1 models because they have different expressions 
for the oscillation probabilities. The difficulty with these 
models is that the sterile neutrino is strongly coupled 
to solar and/or atmospheric neutrino oscillations, and 
there are bounds on the amount of sterile content in each 
case. In 2 -I- 2 models, solar oscillate predominantly 
to a linear combination of Vt and and atmospheric 
oscillate to the orthogonal combination : 

— sin a + cos a J^s , (11) 
cos a Vr + Sinai's . (12) 

If CPT is violated, then there is a similar sterile mixing 
angle a for antineutrinos. The amount of sterile content 
is cos^ a in solar neutrino oscillations, cos'^ a in Kam- 
LAND, sin^ a in atmospheric neutrino oscillations, and 
sin^ a in atmospheric antineutrino oscillations. Fits to 
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solar neutrino data give the 99% C. L. limit 0, 



cos^ a < 0.45 ; 



(13) 



note that there is no limit on cos^ a from KamLAND 
since the short baseline has negligible matter effects and 
it therefore does not test the the sterile content. Fits to 
the atmospheric neutrino data (which do not distinguish 
between neutrinos and antineutrinos) give the 99% C. L. 
limit 0, 



■ sin^ a - 



1 



■ sin^ a < 0.35 . 



(14) 



due to the lack of matter effects that would occur in 
Vs oscillations 26] . In Eq. (|14|l , neutrinos contribute with 
twice the strength of antineutrinos because of their larger 
interaction cross section. 

The bounds on 2 + 2 models are shown in Fig. |2| versus 
sin^ a and sin^ a. The CPT-conserving case a = a is in- 
dicated by the dotted line. A recent comprehensive anal- 
ysis of 2 -I- 2 models with CPT conservation concludes 
that the goodness of fit to all data is only 1.6 x 10~^ 0, 
which is worse than that of 3 -I- 1 models with CPT con- 
servation^. When CPT is violated, i.e., a a, there 
is no region that obeys both the solar and atmospheric 
bounds. Thus even when CPT is violated, 2-1-2 models 
are strongly disfavored. 

The analysis that lead to Eq. (|13|l used the standard 
solar model (SSM) j23| neutrino fluxes, including their 
theoretical uncertainties. If the neutrino flux is al- 
lowed to be free, this bound relaxes to cos^ a < 0.61 [28l| . 
Then the solar constraint in Fig. O is sin^ a > 0.39, and 
there is a small allowed region with 0.39 < sin^ a < 0.53 
and small sin^ a. This can be understood qualitatively 
as follows: although there is no room for a full sterile 
neutrino in the solar and atmospheric data when CPT is 
conserved, with CPT violation the fact that KamLAND 
does not test the sterile content means that oscillations of 
atmospheric antineutrinos can be largely to active neutri- 
nos, which effectively dilutes the constraint on the sterile 
content in atmospheric neutrino oscillations. However, 
there is currently no reason to believe that the neu- 
trino flux is not well-described by the SSM. 



C. Hybrid models 

Since CPT is violated, in principle one can have a 3 -t- 1 
model in the neutrino sector and a 2 -I- 2 model in the 
antineutrino sector, or vice versa. We now examine these 
two possibilities: 



The strong exclusion of the CPT conserving case is also evident 
by adding together the constraints of Eqs. I13i and 1141 . which 
for a = a yields cos^ a + sin^ a < 0.80. 
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FIG. 3: Constraints on sterile neutrino mixing angles a and a 
from solar (solid) and atmospheric (dashed) data. The dotted 
line is the prediction if CPT is conserved. 



(i) (3 -I- l),y, (2 -I- 2)p. In this case, oscillations in the 
neutrino sector are predominantly to active neu- 
trinos, and there is no sterile content bound from 
solar neutrino data. For the atmospheric data, this 
situation is equivalent to a ~ in Eq. (|14|l . which 
leads to no bound on a. Observation of a non- 
negligible sterile content in the oscillations of at- 
mospheric antineutrinos, but not atmospheric neu- 
trinos, would provide supporting evidence for this 
model, although such sterile content is not required. 

(u) (2 + 2)^, (3 + l)p. In this case, oscillations in 
the antineutrino sector are predominantly to ac- 
tive species. The bound in Eq. (|13|l remains the 
same (sin^ a > 0.55), and the bound from the at- 
mospheric data can be found by setting a = 0, i.e., 
sin^ a < 0.53. Therefore this hybrid combination, 
like the pure 2 + 2 model, is strongly disfavored. 

Furthermore, it is possible that more general models 
that do not fit the 3-1-1 and/or 2-1-2 structure can 
provide a satisfactory fit to the data, although we do not 
perform such an analysis here. 



IV. SUMMARY 

We have argued that a four-neutrino model with CPT 
violation can provide an explanation of all neutrino os- 
cillation data if the neutrino sector has a 3 -I- 1 structure 
and the antineutrino sector is either 3-1-1 or 2-|-2;a 
2 + 2 structure is not allowed in the neutrino sector. If 
the antineutrino sector is 3 + 1, then Srnj^ must be near 
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1 or 7 eV^, while if it is 2 + 2, the entire range of (5m| 
allowed by LSND, KARMEN and Bugey is possible. 

A detection of oscillations by MiniBooNE 

with parameters above the solid line in Fig. ^ will re- 
quire CPT violation to be too large and 3+1 models 
with CPT violation will be excluded. However, if Mini- 
BooNE has a null oscillation result or finds oscillation 
parameters that lie below the solid line, such models will 
remain viable. Even if the MiniBooNE allowed region 
for neutrinos falls within the LSND+KARMEN allowed 



region for antineutrinos, the Bugey-|-CDHSW bound on 
CPT conserving 3-1-1 models (dot-dashed line in Fig.0 
implies CPT violation of the size of the LSND effect 
(0(10^'^) in terms of oscillation probabilities). 
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